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1 Introduction 



Experimental pj, I2j and theoretical [3], H] studies of energy loss of fast par- 
ticles provide an extremely valuable option of a quantitative assessment of 
properties of dense strongly interacting medium created in ultrarelativistic 
heavy ion collisions. In such high energy nuclear collisions a transformation 
of incident dense partonic fluxes into a final set of inelastically produced 
hadrons proceeds through initial energy release via creation of physical par- 
tonic modes followed by their transformation into a high density hadronic 
system that finally takes form of free-streaming hadrons that are observed 
experimentally. Detailed dynamical description of high energy nuclear colli- 
sions is an extremely difficult task, so theoretical simplifications and/or com- 
plicated Monte Carlo models are needed to describe the results of existing ex- 
periments and predict those of future ones. The main issue addressed in such 
modeling is a relative role of parton phase contribution versus the hadron 
phase one. One of generic features of a parton phase picture describing dy- 
namics of semihard partonic degrees of freedom is a pronounced turbulence- 
like inhomogeneity of energy release in the impact parameter plane on the 
event-by-event basis [5] having its origin in the powerlike nature of the cross 
sections of partonic interactions. Therefore, when considering the parton 
phase contribution to the energy loss of fast particles, understanding the 
effect of these inhomogeneities is very important. 

In the present paper we shall consider this problem in the simplest set- 
ting [6l [7j in which the energy loss is computed, in the quasi-abelian approx- 
imation, as the work done on the incident particle by the field it induces 
in the medium [8]. In this approach the origin of the losses is a nonzero 
imaginary part of dielectric permittivity of the medium. In the case of ran- 
dom inhomogeneous medium there appears a new contribution to this imag- 
inary part. The physical mechanism of energy loss in the randomly inhomo- 
geneous medium leading to an appearance of such inhomogeneity-induced 
imaginary part of dielectric permittivity can be described as a scattering of 
the Weizsaker- Williams modes of the proper field of the incident particle on 
the inhomogeneities in the medium. Depending on the stochastic proper- 
ties of the medium some modes can become resonant resulting in transition 
and Cherenkov radiation, see e.g. [9]. In the light of discussions on possi- 
ble observation of gluon Cherenkov radiation at RHIC, see e.g. [101 EIj of 
special interest is the fact that the inhomogeneity-induced losses are signif- 
icantly amplified in the medium in which the effective permeability exceeds 
the Cherenkov threshold [9]. Let us note that the problem of energy loss in 
the random inhomogeneous medium has deep links [12] with that of wave 
propagation in random media, see e.g. the review [13] . 
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2 Energy loss 



The average energy losses per unit length of ultrarelativistic particle moving 
along the z axis are described by an appropriate generalization [HI [7] of the 
standard formula 



dW 

dz 



^7Q"z(E"(r,t))r=,i 



where 5^ is a color charge, is a color spin normalized at quadratic Casimir 
invariants of SU(3) for fundamental (adjoint) representations Cp{y^^ , Q°'Q°' = 
Cf(^v), for the projectile quark or gluon respectively and (E°(r, t)) is a chro- 
moelectric field created in the medium by the propagating particle averaged 
over the fluctuations of random inhomogeneous medium. 

In the case of inhomogeneous medium the local fluctuations of its prop- 
erties are customarily parametrized by the coordinate-dependent chromodi- 
electric permittivity e"'''{r). For given configuration of e"'^{r) the spectral 
component of the chromoelectric field is determined from the abelianized 
equations of motion for the chromoelectric field: 



AE" - V (VE'^) + iu^ 6{r) E" = ATiicuf 



(2) 



where r) is a spectral component of the external current and we have as- 
sumed a diagonal structure of chromodielectric permittivity, £:"^(r) = e{r) 5°"^. 
In the considered case of uniformly moving ultrarelativistic particle one has 
j"(co',k) = 2TigQ"''z5{ijj — kz). Random inhomogeneities can conveniently be 
parametrized by explicitly identifying the non-random and random contri- 
butions to dielectric permittivity: 



£(r) = 5o(l + e(r)) 



(3) 



where C, (r) is a random contribution to permittivity having zero mean (r)) = 
0. In what follows we shall consider the simplest case of Gaussian ensemble 
so that ^(r) are fully characterized by the binary correlation function 



(4) 



It is well known, see e.g. [9l|T3], that taking into account the random inhomo- 
geneities leads to an equation for the average electric field containing a tensor 
of effective dielectric permittivity eij(u;,k) having nontrivial imaginary parts 
that depend on the statistical properties of fluctuations ^(r): 



(E>,k)) = — j;>,k) (5) 
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where we have introduced a notation w = ^/eoU. Averaging over the stochas- 
tic ensemble of random inhomogeneities introduces an imphcit dependence 
on its characteristics, e.g. on correlation radius Vc and fluctuation magni- 
tude cr, so that one should in fact write k| Tc, a). It is convenient to 
explicitly introduce the transverse and longitudinal components et and ef. 



£ij{w,k\rc,cr) 



kikj 



kk- 

£*(u7,k|rc,a) + ^£'(w,k|rc,(T) (6) 



In terms of the decomposition ^ the expression for the energy losses takes, 
for some given Sijiu, k| Tc, ct), the form 

1 



dW 

dz 



27r2 

.2 



d'^k 



w 



Im 



£^(w7,k|rc, a) 



(7) 



(w^ — fc^) Im 



w'^e^iw, k| r^, a) — k"^ 



w=k. 



Let us now consider the particular case of an exponential binary correla- 
tion function 



that allows an explicit analytical computation of ejj(u;,k) to all orders in a 
in the one- loop approximation corresponding to the regime a'^{wr^ <C 1, see 
Appendix. The corresponding calculation is naturally performed in terms of 
the polarization tensor njj (w,k) related to chromodielectric permittivity by 
the following relation: 

1 



ei,-(w;,k|rc, a) = eo U 



rlliofw, k r., a. 



(9) 



The decomposition of chromodielectric permittivity (E]) leads to the corre- 
sponding decomposition of the polarization tensor. The explicit expressions 
for the transverse and longitudinal components of the polarization tensor 
njj(w, k| Tc, 0") computed in the one-loop approximation read, cf. Eqs. 
and fl2ip in the Appendix: 

'2 6{6 + iw) 



2 2 

aw 



w 



+ 



w + iSy — k"^ 
6'^ + w'^ + k^ 6 



and 



UHw,k\ rr,,a) 



2 2 

aw 



k^ 



1 + 



-— arctan 
k 



2A;2 
ik 
w + i6 



(10) 



S{S + iw) 



2A;2 



'11^ 



6"^ + w'^ + k'^ 5 



fc2 



-— arctan 
k 



ik 



w + i6 
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where 6 = l/vc- Let us stress once again that even if the non-random chro- 
modelectric permittivity Eq does not have a significant imaginary part so that 
the corresponding energy losses determined by (171) are absent, the effective 
chromodielectric permittivity determined by ( fTOfTTI) has a nontrivial imagi- 
nary part and, therefore, there arise specific energy losses directly related to 
the random fluctuations in the medium in which they propagate. 

Let us now estimate the magnitude of the above-described stochastic 
energy loss. In presenting the results it turns out convenient to rewrite the 
expression for the energy loss (171) in the form: 



where p is an energy of the incident quark (gluon) and the dimensionless 
functions ftiprd c", eo) and fi^prd cr,eo) quantify the transverse and longitu- 
dinal energy losses correspondingly. The dependence on p appeared through 
cutting the momentum integration in ([7]) at p at the upper limit. The phys- 
ical meaning of this cutoff is that within the adopted scheme of calculating 
the energy losses we should restrict our consideration to the soft modes in 
the proper field of the incident particle. 

In considering transverse contribution to the energy losses it is important 
to distinguish between the cases of eo < 1 (below the Cherenkov threshold) 
and Eq > 1 (above the Cherenkov threshold). 

In the first case of < 1 the transverse losses are purely stochastic so 
that they vanish at vanishing fluctuation strength a, ftiprd o" — ^ 0,£o) —* 0. 
An interesting feature of the stochastic energy loss in this regime is a de- 
pendence of the relative weight of transverse and longitudinal contributions 
from the value of Eq. At small Eq longitudinal losses dominate over transverse 
ones while for ^ 1 the transverse losses are dominant. This is illustrated 
in Fig. 1 in which we plot the dependence of transverse ft{prd o"? ^o) and lon- 
gitudinal fi(prd cr,Eo) stochastic losses on pre for e = 0.7, where longitudinal 
losses dominate, and e = 0.9, where transverse losses substantially exceed 
the longitudinal ones. 

In the second case of eg > 1 the transverse losses are nonzero even at 
vanishing fluctuation strength a due to Cherenkov radiation. In the 
considered simplest case of constant energy-independent imaginary part of 
Eo the Cherenkov contribution to the transverse losses is very big: it is in fact 
quadratically divergent at the upper limit of energy integration. Therefore to 
consider the stochastic contribution to energy losses in the Cherenkov regime 
one has to subtract the pure Cherenkov losses from the full transverse ones. 
The stochastic contribution ff^"'^^{prd c", £o) to transverse energy loss in this 
case is thus just ft^°''^{prd cr,Eo) = ft{wd o", ^o) - ftiprd cx = 0, eo)- The full 



dW 





dz 
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transverse loss ftiprd cr, Eq), its stochastic component ft^°^^{prc \ <J, Sq) and the 
longitudinal loss fi^prd cr,£o) are plotted in Fig. 2 for e = 1.1. We see that 
the pure Cherenkov contribution ftiprd u = 0,eo) is clearly dominant while 
the stochastic losses ff^^^ijrrd cr, Eq) and fiijrrd cr, Sq) are similar to those for 
£0 = 0.9 plotted in Fig. 1(b). 

The dependence of transverse losses ftiprd o", ^o) and stochastic transverse 
losses ff'^^'^iprd cr, Sq) on Eq is plotted, for a = 0.2 and = 10, in Fig. 3. We 
see that a relative weight of stochastic losses above the Cherenkov threshold 
is diminishing with growing ^q- 

The dependence of total losses /tot(p'"c| cr, Eq) on the fluctuation magnitude 
a is illustrated by plotting the dependence on pr^ for several values of a in 
Fig. 4 for = 0.7. As expected, the losses grow with growing a. 

To get a feeling on the magnitude of the stochastic energy loss in Gev/fm 
we plot it, for a quark with energy up to 20 GeV and parameter values 
as = 0.2, Tc = 0.2 fm, a = 0.2 and Eq = 0.9, in Fig. 5. We see that at 
high energies the transverse stochastic loss can become quite substantial and 
exceed the values discussed in [HI [E] . 

Let us stress that the considered mechanism of stochastic quenching 
presents and additional source of energy loss with respect to mechanisms con- 
sidered earlier in the literature. It is also important to have in mind that the 
above values constitute an order-of-magnitude estimate of stochastic quench- 
ing only since the above-described calculation is not rehab le at pVc 3> 1, also 
because of the strong dependence of the calculated energy loss on the cutoff 
p. Working out a consistent picture matching the above-described soft en- 
ergy loss to the energy loss of hard modes requires developing a microscopic 
approaclifl. This issue is a very important subject for future analysis. An- 
other important task is to perform a quantitative analysis of the properties 
of inhomogeneous gluon medium created in ultrarelativistic nuclear collisions 
at LHC energies and use it as a basis for computing the stochastic energy 
losses of fast quarks and gluons. . 

3 Conclusions 

Stochastic energy loss of fast color particles in the random inhomogeneous 
medium was calculated to all orders in the fluctuation magnitude in the 
one-loop approximation. The dependence of stochastic losses on the bare 
chromodielectric permittivity Eq, particle energy p and fluctuation strength 
a was analyzed. The analysis has shown that stochastic energy loss can be 

^For an interesting attempt of a combined treatment of ionization and coUisional losses 
see [7]. 
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of phenomenological importance. The strong dependence of the resuhs on 
the value of ultraviolet cutoff indicates a necessity of developing a consistent 
microscopic approach allowing a unified treatment of energy losses due to 
soft and hard modes enabling a reliable quantitative calculation of the cor- 
responding stochastic contribution. It is also of interest to consider in more 
details an interrelation between the above- considered mechanism of energy 
loss and effects of generic stochastic acceleration mechanism in random color 
fields considered in [16]. 
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Appendix: Averaged vector field in statisti- 
cally homogeneous medium. 

To compute an expression for the disorder-averaged chromoelectric field de- 
termining the mean energy loss we have to write down a solution of the 
equation ([2D for some given random profile ^ (r) and then average it over the 
Gaussian ensemble of random inhomogeneities with the two-point correla- 
tor (IH]). The first step is writing down an expression for the Green function 
G(w, r| n, ro) obeying the following equatioru: 

AG - V(VG) + w^G = -w^eWG + nS{r - tq) , (13) 

where n is a fixed unit three- vector and tq is a fixed spatial point. In mo- 
mentum space the equation ffT^ reads 

[{w' - k2)5,, + kik,] G,(k) = -w^iiG^^ik) + jl^n^e-''^''^ , (14) 

where 

(eG)(k) = / dVk'ar)G{k')e<'^'-'^> (15) 



(27r)3 ^ 

Let us consider an iterative solution of ffT^ : 

G^(^'k) = (^G°n,e-^'^--^^^G°^(eX.n.e-^^^^ 

+ (^^'^^^.(^"^^.(e^GL-.-"^'")) + • • • ' (16) 



^In this appendix we omit the color indices. In the considered abehan approximation 
the answer corresponds to a given mode of the chromoelectric field. 
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where G'° (w,k) is a free propagator 

5 in 



+ 



After averaging equation (fTUj) over the gaussian ensemble of fluctuations 
{^(r)} with the two-point correlator ([H]) we arrive at the algebraic Dyson 
equation for the averaged Green function {Gi{w ,'k\ Vc, a, n)) : 



G°,n,e-"^r'' + G°,n,,(G), 



(17) 



so that 



(G(w,k| rc,cr, n)) 



GO(w,k) 



n 



-jkrn 



/ - G^{w, k)n(w, k| r^, a) (27r)3 



The averaged Green function f|T8l) allows to compute the averaged chromo- 
electric field for an arbitrary external current j"(r,t). 

In the limit a'^{krc) -C 1 one can keep the one- loop approximation for the 
polarization operator Tn(w, k| Vc, a): 



Uij{w,k\ rc,a) ^ j d^kig{\ ki - k|, r^, ct)G° (w, ki) 



(19) 



which can be computed analytically. Explicit expressions for its transverse 
and longitudinal components n- (w,k) and n' (t(7,k) take the form 



n*(u;,k|rc,a) 



and 



liHwMre.cr) 



2 2 

aw 



w 



+ 



w + idy — k"^ 
6'^ + w'^ + k^ 6 



k^ 



6{6 + iw) 
2A;2 
ik 



2 2 

a w 



1 + 



-— arctan 
k \w + id 



6{6 + iw) 



2F 



6"^ + w'^ + k'^ 6 f ik 
arctan — 

k \W + 10 



(20) 



(21) 



where 5 = l/r^. 
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Fig. 1(a) 




Fig. 1(b) 
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Figure 1: Transverse losses /t(prc| cr, eo) (red, dashed), longitudinal losses 
fiip^d d'yEo) (blue, dotted) and total losses /tot(p'"c| cr, £^o) (brown, solid), cr = 
0.2: (a) for Eq = 0.7; (b) for Eq = 0.9. 
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Fig. 2 




Figure 2: Full transverse losses ftiprd o", ^o) (brown, solid), stochastic contri- 
bution to transverse losses ft^°'^^{prc\ cr, £o) = ftiprd o", ^o) — ftiprd o" = 0, Eq) 
(red, dashed) and longitudinal losses fi{prc\<J,eo) (blue, dotted), a = 0.2, 

£0=1.1. 
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Fig. 3 




Figure 3: Dependence of full transverse losses ftiprd cr, So) (blue, dotted) and 
stochastic contribution to transverse losses f^^'^'^^^prd ex, £o) = ftip^d o',£o) — 
ftiprd 0" = 0,eo) (red, solid) on the bare dielectric permittivity Eq, cr = 0.2, 
pre = 10. 
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Fig. 4 
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Figure 4: Dependence of total losses ftotip^d o',eo) on pr^ for a = 0.1 (blue, 
dotted); cr = 0.2 (brown, dashed); a = 0.3 (black, dashed-dotted); cr = 0.4 
(red, solid) for Eq = 0.7. 
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Fig. 5 
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Figure 5: Energy dependence of quark transverse (red, solid) and longitudinal 
(blue, dotted) stochastic energy losses, Vc = 0.2 fm, a = 0.2, Eq = 0.9. 
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